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ABSTRACT 

O I Prom an astrobiological point of view, special attention has been paid to the probabil- 

I ity of habitable planets in extrasolar systems. The purpose of this study is to constrain 

' a possible range of the mass of a terrestrial planet that can get water. We focus on 

the process of water production through oxidation of the atmospheric hydrogen — the 
^ . nebular gas having been attracted gravitationally — by oxide available at the planetary 

surface. For the water production to work well on a planet, a sufficient amount of 
hydrogen and enough high temperature to melt the planetary surface are needed. We 
have simulated the structure of the atmosphere that connects with the protoplanetary 
nebula for wide ranges of heat flux, opacity, and density of the nebular gas. We have 
found both requirements are fulfilled for an Earth-mass planet for wide ranges of the 
parameters. We have also found the surface temperature of planets of < 0.3Me (Me: 
Earth's mass) is lower than the melting temperature of silicate (~ 1500 K). On the 
other hand, a planet of more than several Me becomes a gas giant planet through 
runaway accretion of the nebular gas. 



Subject headings: astrobiology — Earth — planets and satellites: formation 
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1. INTRODUCTION 

More than 170 extrasolar planets have been detected so far; most of them are Jupiter-hke 
planets (see www.obspm.fr/planets). At present, several projects are in progress to discover terres- 
trial planets (e.g., TPF, Darwin, etc.). Special attention has thus been paid to habitable planets 
in extrasolar systems. A prerequisite for a planet being habitable is considered to be the existence 
of liquid water on it. For a planet to retain liquid water on its surface, it must acquire a sufficient 
amount of water and be located at a suitable distance from its parent star. This paper focuses on 
the former issue and constrains the probability of the existence of habitable planets in extrasolar 
systems. 

The latter issue has been discussed using the term "habitable zone". The habitable zone (HZ) 
is defined as a range of orbital distance from a star within which a planet can retain liquid water 
on its surface. Around a solar-mass main-sequence star, the HZ is located around 1 AU and its 
width is as narrow as ~0.4 AU (Kasting, Whitmire, &: Reynolds 1993). The surface temperature 
of a planet beyond the HZ is below the freezing temperature of H2O, so that the planet is unable 
to keep liquid water continuously. Because of high surface temperature of a planet closer to its 
parent star than the HZ, the concentration of H2O in the upper atmosphere is so high that the 
planet suffers substantial escape of H2O due to incident stellar UV radiation, and ends up losing 
its ocean completely. 

Formation of terrestrial planets in that narrow zone, however, seems to be by no means un- 
likely within the context of the core accretion model for planet formation. Ida & Lin (2004a, b, 
2005), for example, constructed an integrated model for planet formation including the accretion 
and dynamical evolution of planets. Their model reproduced the mass-period distribution of de- 
tected extrasolar planets and also predicted that of unknown planets. In the predicted mass-period 
distribution, the HZ is filled with hypothetical terrestrial planets of various masses (e.g.. Figure 12 
of Ida & Lin 2004a). This means the existence of planets in the HZ is likely; the remaining issue is 
thus how likely a planet acquires a sufficient amount of water. 

In this paper, we consider the nebular origin of water of terrestrial planets. A planet embed- 
ded in a protoplanetary nebula attracts gravitationally the nebular gas to have a hydrogen-rich 
atmosphere (Hayashi, Nakazawa, & Mizuno 1979; Mizuno, Nakazawa, &: Hayashi 1978). The at- 
mospheric hydrogen can be oxidized by some oxide contained in the planet to produce water on 
the planet, which was first proposed by Sasaki (1990). 

The reason why we focus on the process of water production is that it could commonly happen 
on any extrasolar terrestrial planet. Planets are in general formed in hydrogen-rich protoplanetary 
nebulae. Oxides are also available on a terrestrial planet if the C/0 ratio of the system is less 
than unity (Larimer 1975; Larimer & Bartholomay 1979); most main-sequence stars are known to 
have C/0 ratios less than unity (e.g., Reddy et al. 2003; Takeda & Honda 2005). The remaining 
requirements are a sufficient amount of hydrogen and enough high temperature to melt the surface 
of the planet; the molten planetary surface being called the magma ocean. The second condition 
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is needed because if the planetary surface remains solid, the reaction between the atmospheric 
hydrogen and the surface oxide would result just in production of membrane covering the surface, 
not leading to production of a large amount of water. 

In this paper, we investigate properties of the atmosphere of nebular origin. We first clarify 
the conditions for a planet to have a massive hydrogen-rich atmosphere and a magma ocean to 
constrain the range of the mass of a planet that has a sufficient amount of water. The structure 
of the nebular-origin atmosphere was investigated by Hayashi et al. (1979) and Nakazawa et al. 
(1985) for wide ranges of planetary accretion rate, grain opacity, and density of the nebular gas. 
However, they used quite simple forms of opacity and equation of state, both of which have been 
substantially improved. Furthermore they focused only on the Earth (i.e.. Earth-mass planets) and 
had no discussion on the production of water. Although Sasaki (1990) discussed the production of 
water in order to suggest the deep magma ocean on the early Earth, the ranges of the parameters 
considered were so restricted that we are unable to get any systematic understanding of the nebular 
origin of water on terrestrial planets. 

We also constrain the mass of a planet that remains a terrestrial one. Habitable planets may be 
not gas giant but terrestrial ones. If a planet is isolated from planetesimals, it captures a substantial 
amount of the nebular gas to become a gas giant planet (Ikoma, Nakazawa, & Emori 2000). We 
simulate the evolution and accumulation of the planetary atmosphere to obtain the timescale for the 
substantial gas accretion as a function of planet's mass. The timescale is compared to the lifetime 
of the nebular gas of 10^-10'' years (Natta, Grinin, & Mannings 2000) to constrain a possible range 
of the mass of a habitable planet. 

Section 2 describes our numerical method. Section 3 presents properties of the atmosphere 
of nebular origin for various planetary masses. Section 4 shows the timescale for the substantial 
accretion of the nebular gas. Finally we discuss the probability of water production on terrestrial 
planets and constrain the masses of the potentially-habitable planets in section 5. 



We consider a spherically-symmetric hydrostatic atmosphere that connects with the surround- 
ing nebula. The atmospheric structure is determined by the equation of hydrostatic equilibrium 
including the self-gravity of the atmosphere. 



2. 



NUMERICAL METHOD 



2.1. Basic equations 



dP 

dr 



GMrP 




and the equation of mass conservation. 



dM, 
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(2) 



dr 
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where P and p are respectively pressure and density of atmospheric gas, r is distance from the 
planet's center, is mass inside a sphere of radius r, and G is the gravitational constant. 

The thermal structure is determined by energy transport. If the optical thickness that is 
defined by 

R 

updr (3) 

is smaller than 2/3, temperature is given by (Hayashi et al. 1979) 

A A L 1 + 3t/2 
T^ = T^A ■ '—■ ('4) 

in equations (3) and (4), r is the optical thickness, k is the Rosseland mean opacity, R is the outer 
radius of the atmosphere defined below, T and are respectively atmospheric and nebular gas 
temperatures, and (Tsb is the Stefan-Boltzmann constant. Equation (4) is similar to the well-known 
formula for radiative transfer in a plane-parallel gray atmosphere, but includes roughly the effect of 
spherical geometry. Despite of its roughness, we use this convenient formula because quantities of 
our interest (e.g., surface temperature) are insensitive to the structure of the optically-thin layer. 

If T > 2/3, temperature distribution is determined by the smaller of adiabatic and radiative 
temperature gradients (Kippenhahn & Weiger 1994): The former is given by 

dT GMrP fdT\ 

- ' ' ■ (5) 



dr \ dP 

the latter is given by 

^ SnpL 

dr 64wsBr2r3' ^ ^ 

where S is specific entropy and L is energy fiux passing through a sphere of radius r (called 
luminosity). We assume the opacity is the sum of the Rosseland mean opacities of gas molecules 
and dust grains (see below). The luminosity is determined by (Kippenhahn & Weiger 1994) 

- = -4vrr pT-, (7) 

where we have assumed there is no energy generation other than entropy change in the atmosphere. 

Two types of simulation are done in this paper, quasi-static and static simulations. To calculate 
the timescale of the gas accretion in section 4, we perform quasi-static simulations integrating all the 
above equations but equation (4). Inclusion of equation (4) makes the simulation so complicated, 
but yields little change in the results obtained in section 4. To investigate the properties of the 
atmosphere extensively in section 3, we perform static simulations integrating equations (l)-(6), 
because the quasi-static simulations are time-consuming. In the static simulations we assume L is 
spatially constant, instead of solving equation (7); this assumption was shown to be appropriate 
for less massive atmospheres (Ikoma et al. 2000). 
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2.2. Input physics 

The equation of state used here is the nonideal one given by Saumon, Chabrier, & Van Horn 
(1995), interpolated to a composition of X = 0.74, Y = 0.24, and Z = 0.02, where X, Y, and Z 
are the mass fractions of hydrogen, hehum, and elements heavier than helium, respectively. 

The grain opacity is taken from Pollack, McKay, & Christofferson (1985), who presumed dust 
grains with an interstellar size distribution. Since the amount and size distribution of dust grains 
in the atmosphere are highly uncertain, we regard the grain opacity as a parameter using the 
conventional form of Kgr = /k^^, where the grain opacity, that given by Pollack et al. 

(1985), and / is the grain depletion factor. 

The gas opacity is taken mainly from Alexander & Ferguson (1994), supplemented with the 
data of Mizuno (1980); compared to the latter, the opacity for H2O has been substantially improved 
in the former. In Alexander &; Ferguson (1994), data are available for limited ranges of temperature 
and density, 2.8 < logT < 4.1 and —2 < logQ ^ 3, where Q = p/Tq, p is density in gcm~^, and Tq 
is temperature in millions of degrees. In the calculations here, we need data for lower temperature 
and higher density. For logT > 3.5 both data sets are similar to each other, because a major source 
of the opacity is hydrogen for the temperature range. We thus use Mizuno (1980)'s opacity for 
logQ > 3 and logT > 3.5. For 2.8 < logT < 3.5, a major source of the opacity is H2O, and the 
opacity is rather insensitive to the density. We thus extrapolate values of Alexander & Ferguson 
(1994) to logQ > 3 for that temperature range. For logT < 2.8, we adopt values at logT = 2.8, 
because the opacity is nearly constant for 2.8 < logT < 3.5. 

As mentioned in Introduction, previous workers used rather simple forms of opacity and equa- 
tion of state. Hydrogen molecules dissociate to hydrogen atoms for the ranges of temperature and 
density considered here. Hayashi et al. (1979), however, used the equation of state for a perfect gas, 
although the effect of the dissociation was included in the calculation of the adiabatic temperature 
gradient. The form of the grain opacity adopted by Hayashi et al. (1979) and Ikoma et al. (2000) 
was quite simple; the value of the grain opacity was constant until dust grains evaporate at 1500 K 
independent of gas density. However, the Rosseland mean values of the wavelength-dependent 
opacity of dust grains depend on temperature and, in reality, dust grains consist of several different 
components whose evaporation temperatures depend on gas density. Furthermore, as mentioned 
above, the gas opacity of H2O has been improved and is rather high compared to that used by 
those workers. For low temperature logT < 3.5, H2O makes a dominant contribution to the gas 
opacity. 
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2.3. Boundary conditions 

Four boundary conditions are needed. Two of the four are given at the bottom of the atmo- 
sphere as 

A-Kr^ ( 3Ms \ ^^'^ 

Mr = -^—Ps and L = at r = Rs= [ - — ^ , (8) 
3 V4vrps/ 

where p^-, Mg, and Rg are the mean density (= 3.9 gcm~^), mass, and radius of the sohd part of 

the planet (caUed simply the solid planet hereafter), respectively, and Lg is the energy flux at the 

surface of the solid planet. 

The other two boundary conditions are given at the outer edge of the atmosphere. The outer 
radius is assumed to be the smaller of the Bondi and the Hill radii, which are respectively defined 
by 



1\ GMp 
7/ 4 



Rb={1--]^ (9) 



and 



where Mp is planetary total mass, 7 and ct are respectively adiabatic exponent and isothermal 
sound speed of the nebular gas, a is heliocentric distance, and M^, is mass of the parent star. At 
the Bondi radius the enthalpy of gas is equal to the potential energy of gas, while at the Hill radius 
planetary gravity is equal to the tidal force due to the parent star. For an Earth-mass planet, i?B 
and i?H are approximately 20i?s and 200i?s. 

At the outer boundary the atmosphere is assumed to connect smoothly with the nebular gas, 
namely, 

r = rn and P = Pn at r = R, (11) 

where and are the temperature and pressure of the nebular gas. In this paper, M^, is the 
solar mass, a = lAU, and = 280K, which is the temperature at 1 AU in the minimum-mass 
solar nebula (Hayashi 1981). The nebular pressure (equivalently the nebular density) is regarded 
as a parameter. 



2.4. Parameters 

The parameters in our atmospheric model for a given planet's mass (Mg) are luminosity {Lg), 
grain depletion factor (/), and the nebular density (pn)- 

Because luminosity is supplied by incoming planetesimals in accretion stages, Ls is approxi- 
mately expressed by 

2/3 / M \ / n \ 1/3 



Ls ~ -Ms = 1.1 X 10^'^ [ -] 5 — ^ ^ ^ 5- ergs"S 

Rs \1MeJ I 10-8MEyr-i i V3.9gcm-3; ^ ' 



(12) 
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where Mg is planetesimal accretion rate. The planetesimal accretion rate in late stages of planet 
formation is uncertain. However, in this paper, since we consider the situations where formation of 
terrestrial planets is completed before the nebular gas disappears (10^-10'' years), Mg/Mg should 
be longer than ~ 1 x 10^ years at the time when the planet has grown in mass close to its current 
mass; otherwise, the planet appreciably grows after the disappearance of the nebula. If Mg/Mg > 
1 X 10^ years, Lg < 1 x lO^^ergs^^ for an Earth-mass planet and Lg < 1 x lO^^ergs"^ for a Mars- 
mass (~ 0.1 -Me) planet. A lower limit of luminosity might be constrained by radiogenic luminosity. 
On the early Earth and Mars, the values of radiogenic luminosity are ^ 1 X lO^^ergs-^ (O'NeiU k 
Palme 1998) and ~ 1 x lO^Oergs"^ (Wanke & Dreibus 1988). 

The appropriate range of / is also uncertain. Although Podolak (2003) suggested / ~ 0.01 
through the numerical simulations of coagulation and sedimentation of dust grains in the atmo- 
sphere, more extensive study is needed because their simulations were performed only in a few 
atmospheric models. We thus consider / = 0-1 in this paper. 

The nebular gas dissipates with time. Taking it into account, we consider a wide range of the 
nebular density, 1-10"^" times the gas density (/>msn = 1.2x 10~^gcm~^ at 1 AU) of minimum-mass 
solar nebula model (Hayashi 1981) for the Earth-mass case. 

3. PROPERTIES OF THE ATMOSPHERE 

Figure la illustrates that the surface temperature of an Earth-mass planet changes only by 
at most 50 % despite of the order-of-magnitude differences in luminosity and grain opacity, and is 
always higher than 2000 K when = pmsn- The shapes of the functions for / = and 0.01 are 
convex. The atmosphere is largely radiative for lower luminosities, whereas it is largely convective 
for higher luminosities. Although surface temperature in general increases with luminosity for 
radiative atmospheres, its dependence is found to be weak. This is because the effect of higher 
luminosity is almost compensated for with that of smaller optical depth (i.e., smaller mass of the 
atmosphere) that is due to high luminosity as shown below (see Fig. 3a). For fully convective 
atmospheres, any atmospheric property is independent of luminosity. As luminosity increases, 
the inner convective region expands outward, but, at the same time, the optically-thin, almost 
isothermal layer also extends inward because of a decrease in atmospheric mass. The latter effect 
is responsible for the decrease in surface temperature for high luminosities. The complicated form 
of the function for / = 1 comes from the complicated form of the grain opacity; several different 
components evaporate at different temperatures depending on gas density. 

The surface temperature for / = 1 is clearly lower than that for / < 0.01. To understand the 
difference, we compare the atmospheric structure for / = 1, 0.01, and 1 x 10^'^ in Fig. 2. This 
figure illustrates that the low surface temperature for / = 1 comes mainly from sudden decreases 
in temperature gradient due to evaporation of silicate around 2000 K. For / < 0.01, pressure 
increases substantially in the optically-thin layer and is so high at the evaporation temperature 
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that convection takes place at that point; the evaporation of silicate thus has no influence on the 
structure. 

The values of surface temperature obtained here (< 3300 K) are lower than ~ 4000 K obtained 
by Hayashi et al. (1979). Their analytical calculation for a polytropic atmosphere also yielded ~ 
4000 K. In Fig. 2 we show the atmospheric structure calculated in the similar way to Hayashi 
et al. (1979), using Mizuno (1980)'s gas opacity and the grain opacity of 1 x 10~^cm^g~^ for 
T < 1500 K. The grain opacity of 1 x 10~^cm^g~^ is lower than gas opacity even in the low- 
temperature region. As described in section 2, our gas opacity is much higher than that used by 
Hayashi et al. (1979). Because of the lower opacity in Hayashi et al. (1979)'s model, the isothermal 
layer is deeper, resulting in a large increase in pressure. As a result, in Hayashi et al. (1979)'s 
calculation, the atmosphere is almost fully convective except for the isothermal layer. As shown 
in Appendix, surface temperature for a convective atmosphere is slightly higher than that for a 
radiative atmosphere (e.g., eq. [A17]). This may be the reason why the surface temperature for 
/ = 1 X 10~^ and 1 x 10~^ is slightly lower than that obtained by Hayashi et al. (1979). 

As shown in Fig. lb, the surface temperature of an Earth-mass planet is rather insensitive 
to the nebular density. A decrease in the nebular density certainly lowers surface temperature 
because of a decrease in optical depth, but it changes only by a factor of less than two even if the 
nebular density decreases by ten orders of magnitude. Figure lb shows the surface temperature of 
an Earth-mass planet is always higher than 1500 K, the typical melting temperature of silicate, for 
a wide range of the nebular density. The insensitivity of the surface temperature of an Earth-mass 
planet to luminosity, opacity, and the nebular density is analytically explained in Appendix. As also 
shown in Appendix, it is a rather robust conclusion that the surface temperature of an Earth-mass 
planet is higher than the melting temperature of silicate. 

As shown in Fig. 3a, atmospheric mass decreases almost linearly with luminosity and grain 
depletion factor (for / > 0.01). Lower luminosity and opacity yield gentler temperature gradient 
(see eq. [6]). Then, density gradient must be steeper to maintain sufficiently large pressure gradient 
that supports the gravity. Since the density is fixed at the outer edge, the steep density profile 
results in a massive atmosphere. It should be noted that in Fig. 3a a curve ends off at a critical low 
value of luminosity at which the atmospheric mass is comparable with the planetary mass. Beyond 
the point no static solution is found. That means substantial gas accretion takes place if luminosity 
is smaller than the critical value (see section 4). 

Atmospheric mass decreases as the nebular density decreases, which is shown in Fig. 3b. 
However the dependence is found to be rather weak. Most of the atmospheric mass is concentrated 
in the deep atmosphere. And, as shown in Appendix, the structure of the deep atmosphere is 
insensitive to the outer boundary conditions. 

The atmosphere on a Mars-mass O.IMe) planet exhibits different behaviour compared to 
that on an Earth-mass planet shown above. Figure 4 shows surface temperature (a) and atmospheric 
mass (b) as functions of luminosity for several choices of grain depletion factor (/) and the nebular 
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density {pn)- Unlike the case of an Earth- mass planet, the surface temperature of an Mars-mass 
planet is sensitive to luminosity and grain opacity. The dependence of the atmospheric mass 
of an Mars- mass planet on those parameters is different from that of an Earth- mass planet. In 
Fig. 4a surface temperature is shown to appreciably decrease as luminosity decreases. Although 
atmospheric mass increases with decreasing luminosity in the similar way to the Earth-mass case 
for high luminosities, there is upper limits to atmospheric mass. Comparing Figs. 4a and 4b, we 
find that when the atmospheric mass reaches the upper limit, surface temperature reaches the 
nebular temperature (= 280 K), which means the atmosphere is almost isothermal as a whole. As 
shown in Fig. 4, the surface temperature and atmospheric mass of a Mars-mass planet are also 
sensitive to the nebular density, which is also different from the Earth-mass case. The atmospheric 
mass decreases almost linearly with decreasing nebular density. A mathematical explanation for 
the difference between Earth-mass and Mars-mass cases is given in Appendix. 

Figure 5 shows surface temperature as a function of luminosity for four choices of planet's mass 
(0.3, 0.5, 0.8, and 1 Me), two values of / (0.01 and 1), and = pmsn- Smaller planetary mass 
yields lower surface temperature for given Lg and / (e.g., see eq. [A17]). The surface temperature 
for Ms = 0.8Me exhibits the "Earth-type" behaviour that surface temperature is not so sensitive 
to a decrease in luminosity, while that for Mg = 0.3Me exhibits the "Mars-type" behaviour that 
surface temperature is sensitive to a decrease in luminosity. The case of a 0.5Me planet seems 
to be marginal. As regards melting of the surface of the solid planet that occurs when surface 
temperature is higher than typically 1500 K, the melting is very unlikely on a O.SMe planet even if 
the planet is embedded in the nebula of density as high as that of the minimum-mass solar nebula, 
as shown in Fig. 5. 



4. TIMESCALE FOR GAS ACCRETION 

As mentioned in section 3, no static solution is found below a critical value of the luminosity. 
That means energy supply due to atmospheric contraction is needed to keep the atmosphere in 
hydrostatic equilibrium. The contraction results in substantial accretion of the nebular gas to 
make the planet a gas-giant one. In reality, this happens once the planet becomes isolated from 
planetesimals. In order to exclude gas giant planets from the group of habitable planets, we 
investigate the timescale for the gas accretion (i.e., the timescale for a planet to grow up to be a 
gas giant planet). The gas accretion is known to occur not always in a runaway fashion: Its typical 
timescale increases considerably, as the mass of the solid planet decreases (Ikoma et al. 2000). 

The timescale for the gas accretion also depends on the opacity. As described in section 2, 
Ikoma et al. (2000) used a rather simple form of the grain opacity adopted by Mizuno (1980), while 
we use more complicated and realistic grain opacity given by Pollack et al. (1985). The values of 
Pollack et al. (1985) 's grain opacity are larger by a factor of approximately 10 than those of Mizuno 
(1980)'s for a temperature range from 200 K to ~ 1000 K. Also, unlike Mizuno (1980)'s grain 
opacity. Pollack et al. (1985)'s depends on gas density because the evaporation temperature of rock 
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becomes high with gas density. 

Figure 6 shows the typical timescale for the gas accretion (tg) after planetesimal accretion is 
suddenly terminated; on this timescale a planet becomes a Jupiter-like planet. Following Ikoma 
et al. (2000), we define tg as follows. We first calculate the purely-hydrostatic structure of the 
atmosphere for given values of Mg and Lg (see eq. [8]). Then, setting Lg to be zero, we start quasi- 
static simulation of the evolution and accumulation of the atmosphere. The luminosity emitted 
from the outer edge of the atmosphere soon decreases to a minimum value and, then, increases 
progressively (see Fig. 2b of Ikoma et al. 2000). Most of the gas-accretion phase is spent when 
the luminosity is near the minimum value. The minimum luminosity is almost equal to the above- 
mentioned critical luminosity (Ikoma et al. 2000) . The characteristic time for growth of atmospheric 
mass — Ma/Ma where Ma is atmospheric mass and Ma is gas accretion rate — at the minimum 
luminosity obtained by the quasi-static calculations are plotted with open circles in Fig. 6. Although 
we are unable to obtain the timescale for the gas accretion for / < 0.01 because the simulations for 
/ < 0.01 is quite time-consuming, the result is expected to be similar to that for / = 0.01, because 
difference in / (< 0.01) yields only small difference in atmospheric mass (see Fig. 3a). 

The timescale tg is approximated by (Ikoma et al. 2000) 

GMsM* , , 

tg = a ^ / , (13) 

^conv 

where L* is the critical luminosity, M* and -Rconv the atmospheric mass and outer radius of the 
inner convective layer at L = L*. The curves in Fig. 6 are drawn using equation (13) with the value 
of each quantity obtained by our static calculation and a = 1/3. Those curves are fitted roughly 
by 

for / > 0.01. Not only the value of tg itself but also its dependence on the mass of the solid planet 
obtained here are different from those by Ikoma et al. (2000) in which tg ~ 1 x 10®/(Ms/Me)~^'^ yr. 
In particular tg for Mg = IMe is 100 times as large as that given by Ikoma et al. (2000). This is 
because our grain opacity is larger than theirs. And this is also because our grain opacity depends 
on gas density unlike theirs. For smaller planetary mass, the critical luminosity is small. Then the 
density for a given temperature is high relative to high-luminosity cases, so that the evaporation 
temperatures are also high. Thus the effective grain opacity is higher for smaller planetary masses. 

Similar calculations for a few values of planet's mass were also done by Hubickyj, Bodenheimer, 
& Lissauer (2005). The equation of state and the gas and grain opacity tables they used were almost 
the same as those we have used here. We compare our results with their results of the duration 
(t™'^) between the end of Phase 1 (at which planetesimal accretion was suddenly terminated) and 
the crossover point (at which envelope mass is equal to core mass) in the models named 10H5 
and lOHlO in which core masses are approximately 5Me and IOMe, respectively, and / = 1. 
Although we are unable to exactly compare our results to theirs because the definition of t™^ is 
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not exactly equal to that of tg, both values are found to be shuilar: For Mg = 5Me, tg — 62Myr 
while i™^ ~ 78Myr; for Ms = IOMe, tg ~ 4Myr while t™^ ~ 3Myr. 

5. DISCUSSION 

Based on the numerical results obtained in sections 3 and 4, we discuss production of water from 
the hydrogen-rich atmosphere on a terrestrial planet, a possible range of the mass of a habitable 
planet, and the possibility of the nebular origin of water on the Earth. 

As described in Introduction, production of water on a terrestrial planet requires a sufficient 
amount of hydrogen and surface temperature higher than the melting temperature of silicate (~ 
1500 K). The two conditions are found to be fulfilled on an Earth-mass planet. 

In a late stage of terrestrial planet formation, accretion rate of planetesimals (i.e., luminosity) 
decreases with time because of exhaustion of the planet's feeding zone. The decrease in lumi- 
nosity increases the amount of hydrogen on an Earth-mass planet (see Fig. 3a), while surface 
temperature remains above 2000 K (see Fig. la). For example, if accretion rate of planetesimals 
is 1 X 10~^-/VfEyr~^, corresponding to L ~ 1 x lO^^ergs"^ (see eq. [12]), the mass of atmospheric 
hydrogen is more than about 1 x 10^^ g for / < 1, as shown Fig. 3a: The amount of hydrogen a 
planet acquires is insensitive to the nebular density (see Fig. 3b). 

Water is produced through reaction between the atmospheric hydrogen and oxides contained 
in the solid planet. The amount of water depends on what kind of oxide is available. Ion oxides 
(e.g., wiistite [Feo.9740], magnetite [Fe304], etc.) and fayalite (Fe2Si04) react with the atmospheric 
hydrogen to produce water comparable in mass to hydrogen; the ratios of the partial pressures, 
PH20/PH2, are 0.88, 24.02, and 0.49 at 1500 K for the iron-wiistite (1.894Fe + 02^ 2Feo.9470), 
the wiistite-magnetite (6.696Feo.974 + O2 ^ 2.174Fe304), and the quartz-iron-fayalite oxygen 
(2Fe -|- Si02 + O2 Fe2Si04) buffers, respectively (Robie, Hemingway, & Fisher 1978). Thus, 
if a planet acquires hydrogen of ~ 1 x lO^^g and such oxygen buffers are available, the planet 
obtains water comparable in mass to the current sea water on the Earth (= 1.4 x 10^^ g). However, 
for the silicon-periclase-forsterite buffer (2MgO + Si -|- O2 Mg2Si04), P1120/PH2 is as small as 
~ 3 X 10~^ (Robie et al. 1978). Fe-bearing minerals might be required to produce sufficient water, 
although how much water is needed for a planet being habitable is quite uncertain. Whether Fe- 
bearing minerals commonly exist in extrasolar systems is still a matter of controversy. Equilibrium 
condensation in a highly-reduced environment like a protoplanetary nebula yields not Fe-bearing 
minerals but Fe-metal (Wood & Hashimoto 1993). However, dust grains in a protoplanetary nebula 
can be considered to have non-equilibrium composition including, at least, fayalite (Pollack et al. 
1994, and references therein). 

When the surrounding nebular gas disappears almost completely, the atmosphere and solid 
planet begin to get cold, and then an ocean forms through the condensation of steam in the atmo- 
sphere. However, almost complete dissipation of the nebular gas allows the extremely ultraviolet 
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(EUV) and far-UV radiation from the parent star to penetrate the planetary atmosphere. Such 
irradiation causes extensive loss of hydrogen (and steam). The timescale for complete loss of a 
(say) hydrogen-rich atmosphere due to EUV and far-UV can be estimated to be longer 
than 10^ years (Sekiya, Nakazawa, & Hayashi 1980; Sekiya, Hayashi, & Nakazawa 1981), if very 
strong EUV and far-UV (up to 100 times higher than the present; Guinan & Ribas 2002) from 
the young parent star is considered. On the other hand, the timescale for the ocean formation is 
on the order of 10^ years on a planet located in the HZ, according to calculations based on the 
radiative-convective equilibrium model of an H2O-CO2 atmosphere with mass of ~ 10^^ g (Abe 
1993). The timescale for the ocean formation for a H2O-H2 atmosphere considered here is probably 
not so different from that for an H2O-CO2 atmosphere, because inclusion of H2 hardly affects the 
atmospheric structure due to its weak blanketing effect. Therefore, an ocean can forms before the 
significant loss of the atmosphere due to EUV and far-UV. 

We can constrain a possible range of the mass of a habitable planet. As shown in section 3, 
there is a lower limit to the planetary mass below which the water production proposed in this 
paper does not work. Figure 5 illustrates that surface temperature of a planet of < 0.3 Me is lower 
than the melting temperature of silicate (~ 1500 K) for reasonable ranges of the parameters. On 
the other hand, an upper limit to the mass of a habitable planet can be constrained because a 
massive planet captures a huge amount of the nebular gas to be a gas giant planet, not a terrestrial 
planet. The timescale for the gas accretion depends strongly on the planetary mass, as shown 
Fig. 6. This timescale should be compared to the lifetime of the nebula that is known to be about 
1 X 10'' years (Natta, Grinin, & Mannings 2000). The comparison suggests that the upper limit to 
the planetary mass is 7Me for / = 1 and 2Me for / = 0.01. 

The water production proposed in this paper may have worked on our Earth. Because of 
iV-body simulations of planetary accretion, details of the terrestrial planet formation in the solar 
system have been clarified. After the runaway growth of protoplanets (Wetherill & Stewart 1993), 
they grow in an oligarchic fashion until they eat almost all of the planetesimals in their feeding 
zones (Kokubo & Ida 1998, 2000). Then several Mars-mass protoplanets form in the terrestrial 
planet region. The subsequent growth of the protoplanets needs giant impacts between them. The 
planets formed in the way is likely to have high eccentricities (Ghambers, Wetherill, & Boss 1996). 
Damping of those high eccentricities needs the drag force of the nebular gas (Kominami & Ida 
2002; Nagasawa, Lin, & Thommes 2005). Because the Earth is isolated in the nebular gas, the 
accretion of the nebular gas inevitably takes place and water is produced on the Earth. Although 
the amount of the nebular gas required for the damping of the eccentricities are as small as 10~^ to 
10~^ times that of the minimum-mass solar nebula (Kominami & Ida 2002; Nagasawa et al. 2005), 
our numerical results show that this small amount of the nebular gas is sufficient for the Earth to 
get water comparable in mass with Earth's sea water. 

We should be, however, careful when we consider the nebular origin of water on the Earth. 
This is because the ratio of deuterium to hydrogen (D/H) of the sea water on the present Earth 
is larger by about a factor of seven than D/H of the solar nebula (e.g., Drake & Righter 2002). 
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Moreover, the current Earth's atmosphere includes a tiny amount of noble gas, while the solar 
nebula was rich in noble gas (e.g., Pepin 1991). The latter problem may be solved by extensive 
loss of the noble gases as well as hydrogen from the atmosphere due to EUV and far-UV radiation 
(Sekiya et al. 1980, 1981). Although adequate mixing of water originated from nebular gas and 
water in comets can produce the present D/H on the Earth's ocean, because D/H in comets is 
larger (by about a factor of two) than D/H of the present Earth's ocean (e.g., Drake & Righter 
2002), the former problem is difficult to solve. At present, we have no definite evidences that the 
sea water of the Earth was originated from the nebular gas. 

Our intention in this paper is to claim that water production from the nebular gas on a planet 
is a possible way for a terrestrial planet in the HZ to acquire water. As shown in this paper, the 
water production works, if a planet of 0.3 to several Me forms in the protoplanetary nebula. The 
probability of formation of terrestrial planets in the nebular gas is still open to debate, mainly 
because the dissipation mechanism of the nebular gas is quite uncertain. However, the existence of 
many gas giant planets in extrasolar planets has ensured the validity of the core accretion model. 
That is, accretion of planets generally occurs in a protoplanetary nebula. Also, there is no good 
reason for terrestrial planet formation to prefer vacuum environment. Therefore, it is rather likely 
that terrestrial planets also form in the surrounding nebular gas. The production of water from 
the nebular gas in the HZ thus seems to be a natural consequence of planet formation. 

We are grateful to S. Ida and M. Fujimoto for fruitful discussion and their continuous en- 
couragement. The manuscript benefited from constructive comments from anonymous reviewers. 
We also wish to acknowledge helpful discussion with K. Nakazawa, Y. Abe, T. Tanigawa, and H. 
Senshu. This research was partly supported by the 21st Century COE Program "How to build hab- 
itable planets" , Tokyo Institute of Technology, by Grand-in-Aid for Scientific Research on Priority 
Areas, both of which are sponsored by the Ministry of Education, Culture, Sports, Technology and 
Science (MEXT), Japan, and by the Research Fellowship of the Japan Society for the Promotion 
of Science for Young Scientists. 

A. Analytical solutions of the atmospheric structure 

Analytical consideration of the atmospheric structure gives a deep insight into the atmospheric 
properties shown in section 3. Here we focus on surface temperature and check whether surface 
temperature exceeds the melting temperature of silicate (~ 1500 K). In the following analytical 
consideration, we make some additional assumptions: (1) The atmospheric self-gravity is negligible; 

(2) the atmospheric gas is perfect gas with constant mean molecular weight and adiabatic exponent; 

(3) opacity and luminosity are constant. The first assumption is reasonable because the atmospheric 
mass is negligibly small relative to the mass of a solid planet in cases of our interest. Although the 
second assumption breaks down for T > 2000-3000 K because of dissociation of H2, we adopt it 
because the temperature range of interest is < 1500 K. The third assumption is just for simplicity. 
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but the opacity does not change so much (within less than one order of magnitude) for T < 1500K. 

We adopt the fohowing dimensionless variables, 

P p T r 

w = —, cr = — , = —, and x = —, (Al) 

-TO PO -LQ Kq 

where quantities with suffices are reference quantities. On the above assumptions, using those 
variables, the hydrostatic equation is written by 



1 dw Vq 
a dx 



(A2) 



and equation of state is written by 

w = a9. (A3) 
The equations for convective and radiative energy transports are respectively given by 

w = a^, (A4) 

^ = W,, (A5) 
dw 

In the above set of equations, 

RqPq kRoTo 

where fi is the mean molecular weight, mn is the mass of a hydrogen atom, and k is the Boltzmann 
constant; 

= 4 3 X 10-3 " r - ^ ( ^ ] f^L] fJk_] (A8) 

\IMeJ Vlcm2g-V VlO^^ergs^V VlPay V280Ky ' ^ ^ 

For later discussion, we define an additional quantity A as 

A = ^ (A9) 



7-liis VI^e/ V280Ky V2.34/ V3.9gcm-3y ' V J 



where we have used 7/5 for 7. 

We can obtain the exact solutions for the set of equations (A2)-(A5). The convective solution 
is given by 

while the radiative solution is given by 

l = l + ^{4{e-l)+Fie,Wo)}, (A12) 
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where 



for Wn < 1; 



F{e, Wo) = wo 



In ( e^i±^ V 2 f arctan A _ arctan J- 



Wo 



Wo 



njo^{l-Wo)'/\ 

The exact solution for Wo > 1 can be also obtained (Inaba & Ikoma 2003). 

The dimensionless surface temperature (^s) for a fully convective atmosphere is given by 

1 7-1 



(A13) 
(AM) 



7 



-A, 



since Rq = and thus Vq = 7/(7 — 1); for a fully radiative atmosphere is given by 

e _ ^0 1 _ 7 1 _ A 

' Ax, 4(7-l)xs 4' 



(A15) 



(A16) 



if ^ 1 and <^ 1. The dimensionless surface temperature is found to be determined only by A. 
Thus, surface temperature is given by 



4500 

3900 ( \1MeJ 



IiA/r/ V2.34/ 



Ps 



3.9g cm 



1/3 



K, 



(A17) 



where the upper and lower values are for the convective (7 = 7/5) and radiative cases, respectively. 
Those solutions were also given by Hayashi et al. (1979). The surface temperature given by equa- 
tion (A17) is determined only by planet's mass (for given and ps) and independent of luminosity, 
opacity, and outer boundary conditions. 

The above conclusion is consistent with the numerical results for an Earth-mass planet but in 
contradiction with those for a Mars- mass planet given in section 3. To understand the discrepancy, 
we compare the exact solution (A12) with the approximate one (A16); the result is shown in Fig. 7. 
The difference is found to depend on the value of Wo- For an Earth-mass planet, the difference 
is only 40 % even for logPVo = — 10 and the exact value of Tg ~ 2300 K, while, for a Mars-mass 
planet, the difference is more sensitive to Wq and the exact Tg reaches Tq at logVFo = ~4. From a 
mathematical point of view, the sensitivity is more remarkable for smaller value of A. 

The existence of the outermost isothermal layer lowers surface temperature. To check how 
much surface temperature is reduced, we simulate the structure of a two-layer atmosphere composed 
of the outer isothermal layer and the inner radiative or convective layer. We first determine the 
radius of the photosphere where the optical depth defined by equation (3) is equal to 2/3, using 
the isothermal solution 



a 



exp 



Vn{--1 



(A18) 



Once we obtained values of the quantities at the photosphere, checking convective stability, we cal- 
culate surface temperature using equation (All) or (A12). Figure 8 shows the surface temperature 
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as a function of Pn/ Pmsn for six different values of k and L = 1 x lO^^ergs"^; The solid and dashed 
lines representing convective and radiative lower layers, respectively. If we adopt higher values of 
luminosity, we obtain higher temperatures for k > Ix 10~'^cm^g~^. As shown in Fig. 8, even if 
luminosity is as small as 1 x lO^^ergs"^, surface temperature is above the melting temperature of 
silicate (~ 1500 K) except for a limited range of the nebular density, pn < 1 x 10~^pmsn- 
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Fig. 1. — The surface temperature (i.e., the temperature at the bottom of the atmosphere) of an 
Earth-mass planet for wide ranges of three parameters, luminosity (L), grain depletion factor (/), 
and density of the nebular gas (pn)- In (a)) the surface temperature is shown as a function of L for 
three different values of /. In (b), the surface temperature is shown as a function of normalized 
by that in the minimum-mass solar nebula (pmsn = 1-2 x 10~^gcm~^) for three different sets of 
values of L and /. The attached sets of values are logL in ergs~^ and / in (b). Circles represent 
critical luminosities (see the text). 
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Fig. 2. — The atmospheric structure of an Earth-mass planet for three different values of grain 
depletion factor (/) and luminosity (L) of 1 x 10^^ erg s^^. The structure simulated in the similar 
way to Hayashi et al. (1979) is also shown; in the simulation we have adopted pn = 5.7 x 10~^gcm~^ 
and Tn = 225K, following Hayashi et al. (1979). The bold and dashed lines represent convective 
and radiative layers, respecitively. 
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Fig. 3. — The atmospheric mass of an Earth-mass planet for wide ranges of three parameters, lumi- 
nosity (L), grain depletion factor (/), and density of the nebular gas (pn)- In (a), the atmospheric 
mass is shown as a function of L for three different values of /. In (b), the atmospheric mass is shown 
as a function of pn normalized by that in the minimum-mass solar nebula (pmsn = 1-2 x 10~^g cm~^) 
for three different sets of values of L and /. The attached sets of values are logL in ergs~^ and /. 
Circles represent critical luminosities (see the text). 
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Fig. 4. — The surface temperature (a) and atmospheric mass (b) of a 0.1 Me planet are shown 
as functions of luminosity for several different sets of values of grain depletion factor (/) and 
density of the nebular gas (pn) normalized by that in the minimum-mass solar nebula (pmsn = 
1.2 X 10~^gcm-3). 
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Fig. 5. — Surface temperature as a function of luminosity for several choices of planet's mass. Each 
attached number represents planet's mass in the Earth mass. The solid and dashed lines represent 
cases of / = 1 and 0.01, respectively. 
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Fig. 6. — The typical timescale for substantial accretion of the nebular gas as a function of planet's 
mass for two choices of the grain depletion factor (/). The open circles represent the numerical 
results of our evolutionary calculations; the solid lines are drawn using equation (13) with the 
numerical factor a of 1/3. 
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Fig. 7. — Difference in value of the dimensionless surface temperature (^s) between the exact radia- 
tive solution (eq. [A12]) and the approximate radiative solution called the radiative-zero solution 
(eq. [A16]) for two cases of 1 Me and 0.1 Me planets. The definitions of Wq and A are given by 
equations (A7) and (A9), respectively. 
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Fig. 8. — The surface temperature of an Earth-mass planet as a function of density of the nebular 
gas that is normalized by that of the minimum-mass solar nebula (pmsn = 1-2 x lO^^gcm""^) for 
six different values of opacity (k) and luminosity (Lg) of 1 x lO^^ergs"^ that has been analytically 
obtained in the two-layer (outer isothermal and inner radiative or convective) atmospheric model 
(see Appendix). The solid and dashed lines respectively convective anr radiative inner atmospheres. 



